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Abstract 

The short-time stability of concentration profiles in coherent periodic multi- 
layers consisting of two components with large miscibility gap is investigated 
by analysing stationary solutions of the Cahn-Hilliard diffusion equation. The 
limits of the existence and stability of periodic concentration profiles are dis- 
cussed as a function of the average composition for given multilayer period 
length. The minimal average composition and the corresponding layer thick- 
ness below which artificially prepared layers dissolve at elevated temperatures 
are calculated as a function of the multilayer period length for a special model 
of the composition dependence of the Gibbs free energy. For period lengths 
exceeding a critical value, layered structures can exist as metastable states 
in a certain region of the average composition. The phase composition in 
very thin individual layers, comparable with the interphase boundary width, 
deviates from that of the corresponding bulk phase. 

Keywords: multilayers, Cahn-Hilliard equation, stability, composition, evolution 

1 



I. INTRODUCTION 



Modern techniques of thin film deposition permit the preparation of multilayers with 
nearly arbitrary concentration profiles. The knowledge of the stability of artificial multilayers 
at elevated temperatures is of great practical interest. With increasing atomic mobility, 
compositional changes occur where the diffusion distance is determined by the mobility and 
the available time. In the present paper, the case of coherent periodic multilayers consisting 
of two immiscible components A and B is considered. Accordingly, competing driving forces 
for compositional changes in these multilayers are the reduction of the energy of mixing of 
the two components and the reduction of the interfacial energy of interphase boundaries. 

To evaluate the thermal stability of multilayers, we investigate in the following the 
early stage of compositional changes perpendicular to the individual layers within a one- 
dimensional model ignoring lateral perturbations of the layered structure as well as boundary 
effects at the multilayer surface and the interface to the substrate. The neglect of boundary 
effects is justified as long as the characteristic length of atomic diffusion for the considered 
time scale is small compared to the total multilayer thickness. The compositional evolution 
is studied within the framework of the nonlinear Cahn-Hilliard diffusion equationliHl. Cor- 
respondingly, the multilayer is described by a smooth concentration profile where individual 
layers are separated by diffuse interphase boundaries. The continuum description seems 
questionable in the limiting case of very thin layers of a few monolayers only. However, 
many predictions of the continuum approach agree qualitatively with those of a detailed 
analysis of an appropriate lattice model by HillertS. 

The composition in multilayers evolves quite differently, depending on the initial con- 
centration profile. This is illustrated by the three examples in Fig. |I[ The as-prepared 
multilayers in Fig. [I] were assumed as periodic A/B layer stacks of pure individual layers, 
with the exception of case (c) which exhibits a small thickness perturbation. The curves 
show the concentration of component B as a function of position. They were obtained as 
numerical solution of the Cahn-Hilliard diffusion equation. Depending on the initial layer 
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thicknesses, different cases can be distinguished: (a) dissolution of very thin layers despite 
the immiscibility of the components (phase separation can occur only on a larger length 
scale) , (b) relaxation of a periodic structure to a stationary state with smooth concentration 
profile, and (c) rapid smoothening of the profile followed by slow thinning of thinner layers 
up to their complete dissolution. The main driving force for the composition changes in Fig. 
1 is the minimisation of interface energy. This happens even at the expense of volume energy 
so that the minimal and maximal values of the concentration profile generally differ from 
the concentrations of the corresponding bulk phases. A series of experimental findingsEHil 
reported recently seems to be related to these peculiarities. 

During mechanical alloying, the formation of nonequilibrium supersaturated phases of 
immiscible components as e. g. Ag -Cull and Co-CuBH3 was observed. Due to the small 
layer thicknesses of the lamellar structure which arises during ball milling, a partial inter- 
mixing of the two components could be energetically favourable. Similarly, an enhanced 
solution of carbon in nickel layers was detected in Ni/C multilayers prepared by pulsed laser 
evaporation with individual layer thicknesses of only a few nanometer s0. Another surprising 
observation is the formation of a mixed Co-Cu phase during annealing of Co/Cu multilayers 
despite the immiscibility of Co and The metastable mixed phase formed obviously due 

to the large excess of interface energy in the multilayer. A strong intermixing was also found 
during deposition of a few monolayers Ni onto Aui. In this case, besides interface energy, 
the large elastic energy of the strained Ni layer is an additional driving force for intermixing. 
For a better understanding of these experimental findings, the present theoretical work deals 
with the effect of the high portion of interface energy on the early composition evolution in 
nanoscale multilayers. 

The simulations in Fig. [I] as well as numerical and analytical investigations by other 
authorJE reveal that concentration profiles evolve in characteristic stages: (i) the relax- 
ation to layered quasi-stationary states or the complete vanishing (dissolution) of very thin 
layers takes place comparatively rapidly; (ii) at a longer time scale, a slow ripening process 
involving diffusion between distant layers occurs, i. e. a thinning of the thinnest layers and 



a corresponding thickening of the thicker ones. In the following, the slow ripening process is 
referred to as long-time composition evolution, whereas the relaxation of an arbitrary peri- 
odic concentration profile to a stationary solution of the Cahn-Hilliard diffusion equation is 
referred to as short-time evolution (Fig. [L|b). Also, the rapid dissolution of layers as shown 
in Fig. [TJa is considered as a particular case of the short-time evolution. 

The aim of the present work is to analyse the conditions under which either a rapid 
dissolution of thin layers occurs initially (Fig. [l|a) or a stationary periodic concentration 
profile evolves (Fig. [l|b). To this end, we focus on one-dimensional stationary solutions of 
the Cahn-Hilliard diffusion equation which are characterised by the multilayer period d and 
the average composition c of the multilayer. These two parameters are usually controlled 
by the multilayer preparation and are conserved during relaxation towards quasi-stationary 
states. From the following analysis, a c-d diagram results which shows the existence re- 
gion as well as the stability properties of the stationary solutions. In the present paper, a 
periodic solution is called (globally) stable if its energy is smaller than those of all other 
concentration profiles with the same periodicity d and average composition c. All other 
stationary solutions are metastable or unstable with respect to perturbations conserving the 
periodicity and average composition. As outlined by Langeri, stationary periodic solutions 
are always unstable against perturbations of the periodicity (cf. also Fig. |l]c). However, as 
mentioned above, these periodicity perturbations develop in general on a longer time scale 
and are not considered in this work. Also, small lateral composition perturbations as well as 
a roughening of interfaces between strained inidvidual layers^ are not analysed here because 
it is expected that such perturbations develop on a longer time scale. 

The paper is organised as follows. After the description of the model in Sect. II, station- 
ary solutions of the Cahn-Hilliard diffusion equation are considered in Sect. III. Sect. IV 
deals with a specific model for the Gibbs free energy which allows an analytical calculation 
of stationary concentration profiles. In Sect. V, the stability of these concentration profiles 
is investigated. Finally, the results are discussed and summarised. 
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II. MODEL 



The Gibbs free energy density /(c) of a binary A-B system as a function of the uniform 
concentration c (mole fraction) of component B is to exhibit two minima. The equilibrium 
concentrations, a and (3, of large coexisting phase regions (strictly two half-spaces) are de- 
termined by the common tangent construction (Fig. In the present work, the case of 
multilayers with planar interfaces is considered. Following Cahn and HilliardB, the free en- 
ergy per unit area of a system with concentration c(x) varying in one dimension is described 
by 
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dc 
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(1) 



The second term on the right-hand-side of ([!]) represents the energy contribution due to a 
concentration gradient where k is the gradient energy coefficient. 

The interdiffusion flux in the system is given by j(x) = —M d(5F/ 8c) / dx where M is the 
atomic mobility. Together with the continuity equation dc/dt + Qdj/dx = 0, the following 
nonlinear diffusion equation resultsISi 



- = M — ( SF[C] ' 
dt dx 2 \5c(x) 



M 



dx 2 



f'(c) - 2k 



dx 2 



(2) 



where M = Q M with Q the atomic volume. For simplicity, the atomic volume of the 
two components has been assumed to be equal and the composition dependence of M has 
been omitted. Starting from any initial concentration profile, the further evolution can 
be calculated numerically from @ (see e. g. Fig. 1). However, in view of the great 
practical importance of quasi-stationary concentration profiles, we will analyse the stationary 
solutions of the Cahn-Hilliard diffusion equation (0) more systematically in the following. 



III. STATIONARY SOLUTIONS 



Equilibrium concentration profiles are determined by the extrema of the free energy 
under the constraint of particle conservation. This leads to the variational problem 
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5c(x) 

with the result 



F[c) -fij dxc(x)\ = (3) 



/'(c) -2,^-^ = (4) 

(the prime denotes the derivative with respect to c). Comparison of (f|) and (0) reveals that 
solutions of ([|) are also stationary solutions of the diffusion equation (fj). The Lagrange 
multiplier /i is identified as interdiffusion potential. When \x = SF/5c is uniform, the particle 
flux vanishes. Integration of (f|) leads to 

K {jx) =f(c)-vc + K = D(c) (5) 

where K is an integration constant. The last equality in (|5|) defines the function D(c) used in 
the following. In general, the physically relevant solutions of (j5|) are periodic concentration 
profiles c{x) oscillating between a minimal value a and a maximal value h (cf. Figs. 2 
and 3). The extrema a and h are related to the parameters \x and by the conditions 
D(a) = D(b) = which lead to 

fiP) - /(a) ^ f(b)a-f(a)b 

/x= — , fsT = . 6 

b — a b — a 

Further integration of (|5p yields the inverse function of the concentration profile 

x(c) = [ C dcl{c) (7) 



with 1(c) = yn/D(c) = dx/dc. The integration bounds have been chosen in such a way 
that the origin of the x-coordinate is located at an interphase boundary defined by c = c = 
(a + b)/2. Thus, equation (0) represents the concentration profile in half a period from c = a 
at x — —d a /2 to c = b at x = db/2 (Fig. |3|). The individual layer thicknesses d a and db of 
the two phase regions, briefly called phase 'a' and phase 'b', are given by 

/■co rb 

d a = 2 dcl(c), d b = 2 dcl(c). (8) 

J a J co 

Similarly, the multilayer period length d = d a + db and the mean composition c are given by 
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r b 2 r d b/ 2 2 r b 

d = 2 dcl(c), c=- dxc(x) = - dcl(c)c. (9) 

Ja d J-d a /2 d Ja 

For given concentrations a and b, the concentration profile c(x) can be calculated directly 
from (0). However, from the experimental point of view, a and b are not known a priori. 
Usually, components A and B are deposited consecutively with fixed individual layer thick- 
nesses. During the early annealing stage, a smooth concentration profile develops which 
is similar to the stationary profiles derived here. The arising concentrations a and b are 
determined by equations @ where the mean composition c and period length d are given. 
The solution of equations @ for a and b is, however, not always unique. If there is more 
than one solution, one has to compare the free energies of different solutions in order to find 
that with the lowest one. From (|l|) and ([5]), the free energy of one multilayer period results 
as 

F p = 4v^ [ b dc Jd(c) + 

J a 

An important intrinsic length of the present problem is the width of the interphase 
boundary £ defined by 

e=(6 _ 0) ^i_ =(6 _ a) ^z:. 

The second equality follows from @. In the limiting case of spatially extended phases 
(d a , d b 3> £; i. e. a — > a, b — > (3), the interface width becomes 

£=09-a) V^Zfii =W-a)l (12) 

with f = D((a + /3)/2). f characterises the height of the free energy wall of the /(c) curve 
referred to the common tangent (cf. Fig. |2|). The last equality in (|I^) defines the length 
unit /o used in the following. 

For very thin individual layers, comparable with the width of the interphase boundary, 
the concentrations in the middle of the layers, a and b, differ from those of the corresponding 
bulk phases a and (3 (a < a < b < (3, Fig. |||) because the common tangent construction 
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d . 



(10) 



does not apply to thin layers. The concentrations a and b define a secant with the f(c) 
curve as shown in Fig. |2|. An estimate of the difference between the concentrations (3 and b 
is given by 

f3-b 

- = p b exp(-d h /2£ b ) (13) 

p — a 

(see Appendix A), where = \j2n/ f"(/3) and p b is a numerical factor of the order of unity. 
An analogous formula applies to the difference a — a. In the limiting case d b 3> £&, a factor of 
p b = 2 has been calculated for the special composition dependence of the Gibbs free energy 
considered in the following section. The estimate ( |I~3"D clearly reveals that concentrations j3 
and b differ significantly when the layer thickness d b approaches the characteristic length £&. 



IV. SPECIAL CASE: C 4 MODEL 

To simplify the calculation of the concentration profile (0), we consider the case where 
the free energy as a function of concentration can be represented as a polynomial of the 
fourth power 

f{c) = A{c - af{c - (3) 2 + B{c-a) + C . (14) 

The parameters B and C turn out to be unimportant for the composition evolution. The 
characteristic energy unit fo results as /o = A (j3 — a) 4 / 16. The parameters a and (3 in fljjD 
coincide with the equilibrium concentrations of spatially extended phases corresponding to 
the common tangent construction (Fig. |2|). The characteristic length in ( pT3|) is obtained 
as £ fe = ((3 - a) Z /4. 

For the case flU]) , briefly called 'c 4 -model', the inverse concentration profile (0), the 
period length, and the mean composition (§) can be expressed by elliptic integrals 

x(c)=Uy[j[F e (<!>(c),m)-F e (<P{c ),m)] , (15) 
d = 2uj^K{m) , (16) 
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c = aH — Z(4>z,m) (17) 

with oj = 2/J(b\ — a)(b — m = (b — a)(bi — ai)/((6 — ai)(6i — a)), 0z = 
arcsin — a)/(bx — a{), and 0(c) = arcsin y(6 — ai)(c — a)/{(b — a)(c — ai)). K(m) and 
F e (4>,m) are the complete and incomplete elliptic integrals of first kind, and Z(<fr,m) is the 
Jacobi zeta function@. a and 6 are the concentrations in the middle of the individual layers, 
and at and b\ are the further two intersections of the /(c)-curve with the secant shown in 
Fig. H]. a\ < a < b < b\ are also the four zeros of the function D(c) defined in (|5|). Choosing 
a and b, the other two zeros are given by 

a>i = a + (3 — c — w, bi = a + (3 — c + w (18) 



with w = y 2 [co(a + 0) — ctf3] + ab — 3cq and Co = (a + 6)/2. Analytical expressions for the 
stationary solutions of the Cahn-Hilliard diffusion equation in the case of the c 4 -model have 
been given previously in terms of Jacobian Elliptic functions by Tsakalos0 for the symmetric 
case d a = and by Novick- Cohen and Segel for the general case0. 

Fig. f| shows a series of interesting quantities as a function of the mean concentration 
c for fixed period length d, calculated within the c 4 -model. In a certain c-d region, two 
stationary periodic solutions have been found. The corresponding branches in Fig. ^ are 
denoted by #1 and #2, respectively. The free energy F p of the periodic concentration profile 
with smaller concentration variation b — a (solution 2) is higher than that of solution 1 (Fig. 
|]a). For small values of c, the free energy of the homogeneous concentration, = /(c) d, 
is lower than that of the periodic solution 1. The values of the concentrations a and b are 
shown in Fig. |]b. A striking feature of these plots is the existence of a minimal value c^ of 
the mean concentration. At the minimum off , solutions 1 and 2 merge. The corresponding 
layer thickness of phase 'b' is denoted by d^ 1 (Fig. [|c). It correlates roughly with the 
minimal value of d\, in this case. As discussed in the next section, solution 2 does not appear 
for multilayer period lengths d below a critical value. 
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V. STABILITY OF STATIONARY SOLUTIONS 



In the following, let us consider the stability of periodic concentration profiles c(x) which 
are obtained as solutions of equation ([5]). These concentration profiles are also stationary so- 
lutions of the Cahn-Hilliard diffusion equation @. They are stable against any infinitesimal 
perturbation Sc(x) if the second variation of the free energy 

21 



5 2 F 



dx 



f"(c)(6c) 2 + 2K 



d5c 



(19) 



is positive definite. 

At first, the stability of the homogeneous concentration c(x) = c is considered. Without 



any restriction, the perturbation 5c can be represented as a Fourier series. From (19), it is 
evident that the stabilising influence of the gradient term (second term on the right-hand- 



side of (19)) is stronger the shorter the wavelength of the perturbation is. Considering a 
periodic perturbation with period d, i. e. 5c oc sin(27ra;/<i), we find stability 5 2 F > for 
period lengths 



d < d s (c) 



2vn 



2k 



(20) 



f"(c) ' 

d s (c) is equal to the smallest wavelength of spinodal decomposition obtained from a stability 
analysis of diffusion equation (^)Bi. The inverse function of d s (c) exhibits two branches 
which are denoted by cf(ci) and cf(d) (Fig. ||). For the c 4 -model, one obtains 

1/2 

(21) 



_ s a + /3 p-a 

Cb/aW = — ± 



7T~ 



^b/aK-J 2 - 2v /g 

Within the interval cf < c < cf , the homogeneous solution is unstable and spinodal decom- 
position takes place. According to (0), the maximum of —/"(c) yields a minimal period 
length d m i n below which the homogeneous solution is stable for all values of c. For the 
c 4 -model, d min = tc(P — a) / / v2 follows. 

Since within the region cf < c < cf the homogeneous solution is unstable, there must be 
a stable periodic solution. The investigation of the existence region and stability of periodic 
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solutions is more complicated than that of the homogeneous one. In the following, the 
results of a numerical calculation of the free energies (10) belonging to the solutions flI5|) of 
the c 4 -model are summarised (Fig. ^j. It is expected that the qualitative features are the 
same for other double-well potentials f{c). Similar results have been obtained by HillertB 
for a lattice model. 

The behaviour of periodic solutions changes qualitatively in dependence on the multilayer 
period length d. The stability of concentration profiles can be discussed in a convenient way 
by including non-stationary states with nonequilibrium amplitudes h— a in the consideration. 
Although the parameter b — a does not give a complete characterisation of non-stationary 
states, it is an appropriate quantity to illustrate the stability behaviour of periodic concen- 
tration profiles. In Fig. |6], the free energy F of concentration profiles is sketched as a function 
of their amplitude for different mean compositions and for two multilayer period lengths. 
Extrema of F correspond to stationary solutions. In Fig. [], these extrema are marked by 
filled and open circles denoting stable and unstable stationary solutions, respectively. 

For period lengths above d m i n , but below a certain critical value d c , the situation is re- 
presented by the curves 1 and 2 in Fig. |B|a: (i) Within the region cf < c < cf , there is a stable 
periodic solution corresponding to the minimum of F (curve 1). The homogeneous solution 
(6 — a = 0) belongs to a maximum of F (more precisely a saddle point) and is therefore 
unstable, (ii) Outside the region cf < c < cf , there is no stationary periodic solution (curve 
2). The minimum of the free energy F is given by the homogeneous solution which, therefore, 
is stable. At the concentrations c = cf fe , a continuous transition between the homogeneous 
and the periodic stationary solution occurs. This transition will be referred to as second 
order 'phase transition' in the following despite the fact that the resulting 'phases' are only 
stable in the sense discussed in the Introduction. 

For larger period lengths d > d c , the behaviour is more complex as illustrated in Fig. ^b. 
Within the region cf < c < cf, the situation is the same as for d < d c (curve 1). Outside 
certain limiting values c^f b of the mean concentration, there is no stationary periodic solution 
(curve 4). However, within the regions c^ < c < cf and cf < c < c^f , two stationary periodic 
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solutions exist (curves 2 and 3; cf. also Fig. f|). One solution (solution 1) corresponds to 
a minimum of the free energy F and the other one (solution 2) to a maximum. Thus, 
solution 2 is unstable. Choosing solution 2 as initial condition for the concentration profile, 
the numerical solution of the diffusion equation (fj) reveals a rapid change of the profile. 
Depending on the initial fluctuations, it evolves either into the periodic solution 1 or into 
the homogeneous solution. The latter one corresponds also to a minimum of F. Whether 
the periodic solution 1 or the homogeneous one is more stable depends on the corresponding 
values of the free energy F^ and Fh = f(c)d. At a certain concentration c = c^ 6 (d), for 
which F^ = Fh, a first order 'phase transition' between the periodic and the homogeneous 
states takes place (cf. Fig. (5|). Between c s (d) and c T (d) the homogeneous solution is 
metastable (F^ 1 ' < F h , curve 2 in Fig. ||b), whereas between c T (d) and c M (d) the periodic 
solution is metastable (curve 3 in Fig. ||b). 

At c — c s (d), solution 2 disappears by merging into the homogeneous solution. With 
increasing distance of c from c s , the difference in the free energies of solutions 1 and 2 
decreases and at the concentrations c M they coincide (cf. also Fig. ^). This implies that the 
second variation of the free energy functional vanishes (see Appendix B) and, consequently, 
the periodic solutions become marginally stable at c M . Outside the region c* 1 < c < c^f , 
no inhomogeneous stationary solution exists. Thus, c^f(d) and cff(d) represent the minimal 
and maximal mean concentrations for the existence of metastable or globally stable periodic 
structures with given multilayer period length. 

The critical value of the multilayer period d c , which separates the regions of first and sec- 
ond order 'phase transitions' in Fig. |5], can be derived analytically by means of a third-order 
perturbation analysis. The corresponding critical concentrations c c = c s (d c ) are determined 
by the equation (f") 2 + 3 /" /"" = 0, where /", /"' and /"" denote the second to fourth 
derivative of f(c) at c = c c . According to (p0|), the critical period length is then given by 
d c = 2ii^J 2k/ (— f"(c c )). For the c 4 -model, one obtains c c = [a + [3 ± (/3 — a)/v / 5]/2 and 
d c = \J 5/ 2 d m i n . 

Let us recall that the above picture on the stability of periodic solutions was obtained for 
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the special c 4 -dependence of the Gibbs free energy. Further analysis of other dependencies 
f(c) is desirable to confirm the present stability diagram qualitatively and to study quanti- 
tative changes. Numerical simulations of the composition evolution in the present work were 
performed by means of a finite difference method described in ref.0. This method is based 
on a semi-implicit finite difference scheme coupled with a fast Fourier transformation. For 
numerical details, the reader is referred to Copetti and The spatial grid spacing 

in our calculations was less than 10% of the interphase boundary width of stationary states 
and for the time integration an adaptive step size control was applied. For example, 1024 
grid points were used for the calculations in Fig. 1. 



VI. DISCUSSION 

Multilayers for practical applications are usually deposited as layer stacks of nearly pure 
layers of components A and B with thicknesses and d>B- During annealing, a considerable 
interdiffusion between individual layers can occur depending on the values of the equilibrium 
bulk concentrations a and (3, which change with temperature. Consider for definiteness the 
case of thin layers of component B between thick layers of component A (ds < d^). For 
arbitrary multilayer period lengths, the B-layers dissolve at elevated temperatures if the 
mean composition c = ds/d (supposing Qa — Rb) is smaller than a, where a increases with 
increasing temperature. However, layer dissolution can occur also for c > a, if the multilayer 
period is small enough. 

The composition evolution in multilayers is controlled by the competition between volume 
free energy and interface energy. For mean compositions outside the region c^f(d) < c < 
c^(d), no stationary periodic solutions of the Cahn-Hilliard diffusion equation have been 
found. This implies that individual layers of such multilayers dissolve in the early stage of 
annealing because the free energy gain by phase separation is too small to overcompensate 
the interface energy (Fig. la). The thickness of the 'b'-layers d^f = db(c = c* f ), which 
corresponds to the minimal mean concentration c^f(d), is shown in Fig. [7| as a function of 
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the period length d. For large period lengths (d > 7 Iq for the example in Figs. [| and 
d^f(d) represents the minimal layer thickness of phase 'b' of stable concentration profiles. 
For smaller d (but larger than d c ), the minimal thickness dj, appeared in the present case at 
a mean concentration slightly larger than c^ 1 (cf. Fig. |]c). 

Despite the large change of the multilayer period length in Fig. [7|, the characteristic layer 
thickness d^f(d) changes only slightly and is always of the order of the characteristic length 



= y K /fo- The gradient energy coefficient k can be estimated within the framework of the 
regular solution mo deffl AU/r where r is the interatomic distance and AU is the 

energy of mixing (per atom) of an equiatomic solution (c = 0.5); typically, « = 10~ n to 10~ 10 
J/m. The value of Iq considerably exceeds the interatomic distance when / <C AU/tq, where 
AU/tq is typically in the range of 10 8 to 10 9 J/m 3 . Choosing for example k = 3 ■ 10~ n J/m 
and fo = 10 8 J/m 3 , one obtains Iq ~ 0.5 nm which is about twice the interatomic distance. 

The gradient energy coefficient changes only very weakly with temperature, whereas the 
parameter f (cf. Fig. 2) decreases with increasing temperature, approaching zero at the 
critical phase separation temperature T c . As a consequence, the length l and correspond- 
ingly the thickness dff, characterising the onset of layer dissolution, diverge as \T — T c | -1 / 2 
at the critical temperature. Comparatively small critical temperatures in technologically 
important regions of a few 100°C are found, for example, for systems with large lattice mis- 
match because mechanical stresses in coherent layers cause a considerable lowering of the 
critical temperature compared to incoherent phases§l. 

The layer dissolution in multilayers with very small period length, which is driven by the 
reduction of interface energy, leads to the formation of supersaturated phases. Subsequent 
phase separation occurs with a larger period length on a longer time scale and can be kine- 
tically hindered by rapid quenching. As mentioned in the introduction, such a situation could 
be present during mechanical alloying^, when a nanoscale lamellar structure develops in 
the course of ball-milling, or in the case of multilayer deposition0lll with subsequent short- 
time annealing at moderate temperaturesEl. The nonequilibrium phase formation observed 
in those experiments could be related to the layer dissolution discussed here. Based on an 
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estimate of the chemical energy of coherent phase boundaries, such a mechanism has been 
suggested by Gente et al.l to explain the observed solid solution formation of immiscible 
elements due to long-time ball-milling. 

VII. SUMMARY 

Stationary solutions of the one- dimensional Cahn-Hilliard diffusion equation have been 
analysed in dependence on the mean composition c of a binary system and the multilayer 
period length d. A c-d diagram has been established showing the regions of existence, 
metastability and global stability of stationary periodic concentration profiles, as well as of 
the homogeneous concentration (Fig. ||). The diagram was derived under the constraint of 
fixed period length of the concentration profile. Actually, periodic solutions are unstable 
against thickness fluctuations! which leads to layer thickness coarsening in the course of 
annealing. However, this happens on a longer time scale than the relaxation to quasi- 
stationary concentration profiles as long as the individual layer thicknesses are significantly 
larger than the interphase boundary width. 

The present analysis revealed that very thin individual layers (typically a few monolayers 
for small mutual solubility) dissolve during annealing if the mean composition of artificial 
multilayers is lower than a critical value c^(d). The individual layer thickness d^ 1 corre- 
sponding to the minimal value c^(d) increases slightly with increasing multilayer period 
length. The layer dissolution is driven by two mechanisms: (i) interdiffusion between pure 
individual layers to establish the equilibrium phase concentrations and (ii) reduction of in- 
terface energy in the case of very thin layers of the order of the interphase boundary width. 
In a certain c-d region, the layered structure can exist as metastable state although the free 
energy of the homogeneous concentration is lower. According to the Cahn-Hilliard theory, 
the equilibrium composition in thin layers differs from that of the corresponding bulk phase 
if the layer thickness becomes comparable with the interphase boundary widtrJHl. 

Consideration of the present conclusions in the design of layered structures could help 
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to improve their thermal stability or, on the other hand, to prepare new metastable phases 
by controlled layer dissolution. Without changing the qualitative predictions, mechanical 
stresses in coherent layers due to lattice mismatch can easily be included in the present 
one-dimensional analysis by a modification of the free energy f(c) (see e. g. ref.i). 

The evaluation of the long-time stability of multilayers requires an additional analysis of 
the evolution of lateral composition perturbations including the effect of stresses in individual 
layers. The roughening instability of interfaces between strained layers has been investigated 
for example in refJil. A comprehensive analysis of the composition evolution under the 
influence of stresses owing to lattice mismatch or to the presence of dislocations has been 
given recently in a series of papers by Leonard and DesaiS. 
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VIII. APPENDIX A 



To derive (]T3|), equation (f|) is linearised with respect to the concentration difference 
Scb(x) = b — c(x). With the approximation f"(b) ~ f"(P), one obtains 

[d^-eJ^^Hr 1 (22) 

(£ 2 = 2k/ f"(f3)). Choosing the middle of the 'b'-layer as origin of the x-coordinate and 
requiring 5q,(0) = 5c' b (0) = 0, the solution of (|22"D is obtained as 



Although this result was derived for small values of 5cb, it is extrapolated to larger values in 
order to get an estimate for the layer thickness db defined by c(±c4/2) = (a+b)/2. Assuming 
further db 3> £b an d approximating b — a by (3 — a, one finds 

Scb(±db/2) = Lz£ „ Wzfi = I <lz/|) expW/2&) . (24) 



The correction factor is introduced in fl2l|) in order to account for the error caused by 
the linearisation of ([22] ) within the interface. The last equation can be rewritten using 
the expansion f'(b) ~ f'(P) + f"(P)(b — (3) and the fact that \i — f'(/3) is of the order 
of 0((/3 — b) 2 , (a — a) 2 ). Neglecting these higher order terms and inserting \x — f'(b) ~ 
f"(P)(P — b) into (H), one obtains equation dl3|) . 



IX. APPENDIX B 

In the following, the stability of two stationary solutions for the same d and c, c\(x) 
and c 2 (x) differing by 8c(x) = c 2 (x) — ci(x), is analysed for 5c — > 0. The two solutions 
fulfil equation (f|) with the corresponding Lagrange multipliers H\ and /i 2 . Expansion of (|J) 
for C2 = ci + 5c with respect to 5c and 5fi = fi 2 — A*i yields the following equation for 5c: 
f"(ci) 5c — 2Kd 2 5c/dx 2 = 5fi. Using the estimate 5ji ~ 0((5c) 2 ), one obtains to first order 
in 5c 
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/"(ci) - 2k 



dx 2 



Sc(x) = . 



(25) 



This equation is equivalent to marginal stability of c\{x). Indeed, the second variation ( |T9|) 
of F[c] can be transformed by partial integration into 

r T d 2 ' 

5 2 F = dx5c f"{c) - 2k— 5c . (26) 




In deriving (^), the periodicity of solutions C\{x) and 02(2), and consequently of 5c(x), was 

used. Comparison of (|26|) and (|25|) leads to 5 2 F = 0. In summary, if two solutions merge, 
marginal stability results. This happens at the boundaries c s as well as c M (Fig. ||). 

* Present address: Max-Planck-Institut fur Physik komplexer Systeme, Nothnitzer Str. 
38, D-01187 Dresden, Germany 
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FIGURES 

FIG. 1. Early composition evolution of stacks of initially pure layers (t = 0, dashed lines): (a) 
rapid dissolution of individual layers (t/to = 0.004, 0.02, 0.1), (b) relaxation to quasi-stationary 
profile (t/to = F 1000; the two curves coincide), and (c) fast relaxation to quasi-stationary profile 
with subsequent slow thinning and rapid dissolution of the initially thinner layer (t/to = 0.1, 171.3, 
171.5; to = k/M/q, parameters see Sect. II and III). 

FIG. 2. Gibbs free energy density as a function of concentration of component B together with 
the common tangent and a secant corresponding to infinite and finite period length d of stationary 
concentration profiles, respectively. Fhe secant intersects the f(c) curve at the concentrations 
a\ < a < b < b\. 

FIG. 3. Stationary concentration profile illustrating the meaning of concentrations a and b and 
the layer thicknesses d a and db (a = 0.1, (3 = 0.9). The dashed line represents a second stationary 
profile for the same parameters d and c, which is, however, unstable (cf. Sect. V). 

FIG. 4. (a) Free energies of periodic concentration profiles F p (full line) and of the corre- 
sponding homogeneous concentration Fh (dashed), (b) concentrations a (full) and b (dashed), as 
well as (c) layer thickness of phase 'b' as a function of the mean concentration c (d = 6lo, a = 
0.01, j3 = 0.99). The intersection of F p (c) and FX(c) determines the phase boundary c T in Fig. ||. 
The thickness db at the minimal value c^f defines d^ in Fig. 0. 

FIG. 5. Phase diagram showing the phase boundary c T (d) between homogeneous and periodic 
solutions as well as the limits of metastability of the periodic and homogeneous solutions c M (d) and 
c s (d), respectively (a = 0.01, f3 = 0.99). The vertical dashed line at d = d c separates the regions 
of first and second order transition at the concentrations c s (2-PT) and c T (1-PT), respectively. 
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FIG. 6. Schematic of the dependence of the free energy on the amplitude of the concentration 
profile for two period lengths (a: d < d c , b: d > d c ). The origin of the amplitude-axis corresponds 
to the homogeneous concentration b — a = (for further explanations see text). 

FIG. 7. Layer thickness df = <4(cf ) as a function of the period length d; in the unstable 
region individual layers dissolve (a = 0.01, (3 = 0.99). 
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